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ABSTRACT 


This  thesis  investigates  the  multiplexing  of  Mach-Zehnder  type  3x3  terminated  fiber¬ 
optic  sensors  demodulated  by  either  'quadrature'  or  'symmetric'  methods  using  intensity 
modulation  of  the  source.  3x3-couplers  produce  signals  that  permit  unmodulated  passive 
demodulation  of  interferometric  signals.  The  theory  is  described  and  the  results  of  a  2x1 
element  array  optical  demonstration  are  presented.  Possible  architectures  using  this 
demultiplexing  technique  are  presented  for  several  applications  with  different  return  line 
requirements.  The  technique  was  successfully  demonstrated  and  warrants  further 
investigation  to  increase  the  number  of  sensors  and  reduce  the  number  of  return  lines  for 
specific  applications.  The  multiplexing  technique  presents  the  opportunity  for  possible  cost 
savings  over  other  phase  generated  carrier  techniques,  which  require  wavelength 
modulation  of  the  source,  and  significant  optical  path  differences  in  the  interferometers, 
and  are  therefore  constrained  to  presently  very  expensive  sources.  The  technique 
presented  uses  compatible  low  coherent  laser  sources  such  as  Compact  Disc  quality  (830 
nm)  devices. 
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TABLE  OF  SYMBOLS 


Symbol 

Definition 

A 

The  sum  of  the  squares  of  the  two  interfering  electric  fields  and  will  be  termed 
the  dc  bias  value  of  the  optical  intensity. 

B 

Twice  the  product  of  the  two  interfering  electric  fields  which  reprents  the 
amplitude  of  the  cosinusoidally  varying  signal  whose  argument  is  the 
differential  phase.optical  intensity  in  the  fiber,  i.e.,  the  maximum  variation  from 
the  dc  value. 

C 

The  amplitude  of  the  mixing  sinusoidal  signal. 

c„ 

The  coupler  ratio  for  coupler  n  The  coupler  ratio  is  the  power  split  ratio  of  a 
particular  coupler  as  a  fraction  of  the  total  input  power,  i.e.,  a  percentage  or 
equivalent  fraction. 

F 

The  number  of  return  lines  for  a  given  array  architecture. 

G 

Gain  factor. 

Jn 

The  Bessel  function  of  order  n. 

k 

The  wave  number  in  a  vacuum  defined  as  k  =  2nlX,  where  X.  is  the 
wavelength  in  meters  unless  stated  otherwise. 

L 

The  physical  length  of  optical  fiber  in  meters  unless  stated  otherwise. 

Lc 

The  coherence  length  of  a  laser,  defined  as  the  distance  over  which  the 
amplitude  modulation  falls  off  to  an  arbitrarily  determined  fraction  of  the 
complete  modulation  value  [Reference  6,  page  46]. 

M 

The  number  of  lasers  in  an  array  of  sensors. 

m 

The  modulation  index,  a  measure  of  the  amount  of  modulation  as  a  fraction  of 
the  total  available,  having  a  value  between  0  and  1. 
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Symbol 


Definition 


//,  V,  and  w 


The  number  of  sensors  per  laser  in  an  array  of  sensors. 


The  refractive  index  of  a  particular  medium  where  the  refractive  index  is 
defined  as  the  ratio  of  the  speed  of  light  in  a  vacuum  to  the  speed  of  light  in  the 
medium. 


Pfj  denotes  the  output  power  from  a  particular  laser.  Subscripts  are  used  to 
denote  the  power  from  a  particular  coupler  in  an  array,  e  g.,  the  power 
transmitted  by  the  coupler.  P'  denotes  the  power  normalized  to  Pq. 


The  differential  phase  shift  in  the  arms  of  the  interferometer.  In  general  ^  is  a 
time  dependent  quantity,  however  the  time  dependent  notation  is  dropped 
for  compactness.  It  is  often  used  with  various  subscripts  to  denote  the  origin  of 
the  phase  shift  e  g.,  <l>(jrij^t)  denotes  the  phase  shift  due  to  factors  other  than 
the  signal  of  interest.  denotes  the  amplitude  of  the  sinusoidal  signal  of 
interest. 


The  amplitude  of  wavelength  modulation  which  occurs  at  the  modulation 
frequency. 


The  total  number  of  sensors  for  a  given  array. 


The  frequency  of  the  phase  shift  in  the  interferometer,  either  in  Hertz  or 
Radians  respectively. 


The  light  intensity  of  each  of  the  three  legs  of  a  3x3  terminated  Mach-Zehnder 
interferometer.  Each  is  phase  shifted  120°  relative  to  the  other  two. 
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1.  INTRODUCTION 


.1.  MOTIVATION 

The  United  States  Navy  is  seeking  to  take  full  advantage  of  fiber  optic  technology,  in 
particular,  the  bringing  together  of  fiber  optic  conununications  technology  and  fiber  optic 
sensor  technology.  Fiber  optic  communications  systems  have  many  advantages  over  other 
methods  of  data  communication.  These  include  high  bandwidth,  low  power  consumption, 
security,  electromagnetic  interference  (EMI)  protection,  electromagnetic  pulse  (EMP) 
protection,  light  weight,  low  loss,  and  low  cost.  Optical  fiber  sensors  based  on 
interferometry  have  been  designed  with  high  sensitivity  and  dynamic  range  Various 
demodulation  and  multiplexing  schemes  are  feasible  for  use  with  these  interferometric 
sensors  [References  1,  2,  3,  4,  and  5].  For  example  it  has  been  demonstrated  that  3x3- 
coupler  terminated  interferometric  sensors  with  quadrature  [Reference  6]  and  'symmetric' 
demodulation  yield  particular  advantages  [References  7,  8,  and  9].  An  array  of  these 
optical  fiber  sensors  in  an  optical  architecture  has  the  potential  to  accommodate  both  the 
benefits  of  the  interferometric  sensors  and  optical  communications  systems.  An  all-optical 
system  removes  the  need  for  outboard  electronics  and  thus  has  the  potential  for  great 
weight  and  cost  savings.  Such  an  optical  system  using  3x3  couplers  and  'symmetric' 
demodulation  is  demonstrated  in  this  thesis. 

B.  GOAL 

This  thesis  demonstrates  that  the  use  of  3x3 -coupler  terminated  interferometric 
sensors  and  either  quadrature  or  'symmetric'  demodulation  enables  the  use  of  amplitude 
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modulation  (AM)  frequency-division-muftiplexing  (FDM)  techniques  where  the  laser 
intensity  is  the  quantity  that  is  amplitude  modulated 

The  theory  and  performance  of  such  a  multiplexing  scheme  is  investigated  with  sensor 
demonstrations  using  both  analog  electronics  and  optical  fiber  interferometric  sensors  On 
the  basis  of  this  research  a  design  for  an  array  of  M  x  N  sensors  is  presented  which  may 
provide  large  cost  savings  with  improved  system  performance. 


C  BACKGROUND 
1.  Interferometric  Sensors 

Interferometry  is  a  technique  used  to  measure  differential  phase  shifts  The  optical 
phase  delay  generated  in  a  physical  length  of  fiber  L  is. 

<f>  =  hiL,  (11) 


where  n  is  the  refractive  index  of  the  fiber  core  and  k  is  the  wave  number  in  a  vacuum  (k  - 
2x/X).  Neglecting  birefringence  from  the  two  orthogonal  polarizations  supported  by  the 
fiber,  the  fractional  variation  in  phase  can  be  expressed  as; 


d<t> 


k 


dn 


n 


dL 
L  ’ 


(12) 


where  the  changes  in  n  and  L  are  related  to  a  particular  measurand  (e  g.,  pressure, 
temperature,  magnetic  field,  or  other  environmental  changes). 

Various  techniques  have  been  demonstrated  to  selectively  enhance  the  change  in 
phase  in  the  fiber  for  a  particular  measurand  (e  g.,  acoustic  waves),  while  reducing 
changes  in  phase  due  to  unwanted  effects  (eg.,  acceleration).  Mach-Zehnder 
interferometric  sensors  use  two-beam  interferometry  as  illustrated  in  Figure  1.1.  It  has 
been  shown  that  by  using  both  legs  as  signal  legs,  in  push-pull  configurations,  improved 
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sensitivity  and,  through  design,  highly  favorable  properties  (e  g.,  acceleration  canceling  ) 
can  be  obtained. 


Figure  1.1.  Schematic  diagram  of  a  Fiber  Optic  Mach-Zehnder 
Interferometric  Sensor 


An  ideal  3x3-coupler  provides  three  signals  of  equal  amplitude,  each  with  a  relative  phase 
difference  of  120°.  The  intensities  of  each  output  leg  of  the  3  x  3-coupIer  are: 

u  =  A  5cos  +  (13) 


Bcos  [<^{t)  +  (f>,„j,{t)  +  \20°\ 

(1,4) 

Bcos 

(1.5) 

where  A  is  the  value  of  optical  intensity  from  the  fibers  (analogous  to  a  dc  offset),  B  is  the 
peak  value  of  the  intensity  and  is  the  static  phase  shift  resulting  from  the  difference 

in  optical  path  length  between  the  two  signal  legs.  This  quasi-static  (or  slowly  varying 
compared  to  the  signal  of  interest)  phase  also  includes  the  phase  shift  due  to 
environmental  effects  (hydrostatic  pressure,  temperature,  etc.)  as  well  as  the  physical 
difference  in  the  length  of  the  legs.  Examples  of  these  types  of  signals  are  shown  in 
Appendix  A.  The  signal  of  interest  <j)(t)  is  assumed  to  be  of  the  form; 
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^(0  =  ^,sin(<y,0 


(16) 


Using  the  trigonometric  identity: 

cos(a+>3)  =  co^(a)cos(^+sin(a)sinO!?),  (17) 

the  intensity,  equation  (1.3),  can  be  rewritten  as; 

M  =  ^  +  5cos[^,sin(fl),/)]cos[^^^(/)+^,] 

-  fisin  [<!>,  sin  {a,  /)]  sin  [^rf„y,(/)+ 

Similar  expressions  for  v  and  w  with  phase  shifts  of  ±120°,  as  appropriate,  also  can  be 
written.  This  equation  is  of  the  form  of  a  Bessel  function  generating  function.  The 
generating  function  and  associated  series  from  Reference  10  (page  361,  equations  9  1  42 
and  9.1.43)  are  reproduced  below; 

ao 

cos(z  sin  ff)  =  Jq  iz)  +  2^  (z)  cosilkff)  ( 1 .9) 


sin  (r  sin  =  ,(z)sin[(2A:  +  l)4  (1  10) 

4=0 


The  intensity,  equation  (6),  can  therefore  be  written  as  a  series  of  Bessel  weighted  sine  and 
cosine  functions  as  follows. 


u  =  A  +  Bcosi(^j„^it))\ 


lr=l 


(1.11) 


-^sin(^,,„^(/)) 


2^  (^.  )sin  [(2A:  +  l)fl;,/] 

.  *=o 


Thus  the  frequency  spectrum  consists  of  components  at  multiples  of  co^  which  are  in  turn 
weighted  by  Bessel  function  terms  and  the  drift  term,  Bcos{^j„j,(t)).  Appendix  A  plots 
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the  lower  order  Bessel  functions  and  the  frequency  spectrum  of  Equation  1. 1 1  for  typical 

HO  and  H„fi(0 

2.  Symmetric  Demodulation 

The  technique  of  symmetrically  demodulating  the  outputs  was  developed  by 
Professor  Robert  M.  Keolian  of  Naval  Postgraduate  School  and  is  detailed  in  References 
7,  8,  and  9.  A  schematic  of  the  scheme  is  shown  in  Figure  1.2.  Modifications  and 
performance  specifications  of  the  symmetric  analog  demodulator  of  Reference  9  are 
detailed  in  Appendix  B. 


Figure  1.2.  Passive  Symmetric  Demodulation  Scheme  (adapted  from  Reference  8). 


The  demodulation  scheme  is  based  upon  the  following  expression; 


4>(t) 


f  u(v- w)-^  v(w-u)-^w(u~  v) 


J- 


(u^  +v^  +  w^) 


dt. 


(1  12) 


where  ^i)  is  the  signal  of  interest  to  be  recovered  and  the  dots  above  a  variable  indicate  a 
time  derivative.  The  implementation  detailed  in  Reference  9  is  outlined  below  for 
continuity. 

The  output  intensity  of  each  leg  of  the  interferometer  is  as  previously  defined; 

u  =  A  +  5cos[^)(0  +  ^</„/j(0].  (113) 

with  similar  expressions  for  v  and  w.  After  photodetection  and  amplification  by  low-noise 
transimpedance  amplifiers,  the  symmetric  voltages  are  added  together  and  scaled  by  a 
factor  of  -1/3  to  produce  -A.  This  average  dc  content  signal  is  then  summed  to  each  of 
the  symmetric,  interferometric  voltage  signals  to  remove  the  average  dc  component,  A. 
These  interferometric  signals,  with  the  average  dc  component  removed  are  denoted  again 
by  //,  V,  and  w,  where 

u  =  5co^^/)  + 

V  =  5cos[<^(/)+^,„/t)  +  120°]  (1.14) 

W  =  5cOs[$K0+^,fr,y»  (0-120°]. 

These  three  symmetric  interferometric  signals  are  cross  multiplied  by  the  difference 
of  the  'complimentary'  derivatives  to  yield; 
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u{v-w)  =  5cos(^>)|^5^sin(^-120")  +  5^sin(^  +  120°)j 

v{w-u)  =  5cos(^-120°  sin(^  +  120®)  + sin(^)j  (115) 

w(m-v)  =  ficos(^-120“)|^5^sin(^  +  120“)  +  5^sin(^)j. 

It  can  be  shown  by  trigonometric  identities  or  with  phasor  diagrams  that  the  first  equation 
of  (1 . 16)  can  be  simplified  to; 

m(v-h')  =  y/3B^  ^  cos^ (^) .  (1  16) 

By  noting: 

cos‘(x)  +  cos*(jf  +  120*)  +  cos‘(x-120°)  =  4,  (1  17) 

the  sum  of  the  three  equations  in  (1 . 16)  can  be  shown  to  result  in: 

Ji5'[cos'(^)  +  cos'(^+120^)  +  cos2(^-I20®)]  =  (1.18) 

The  factor  of  B^  varies  as  a  function  of  laser  intensity,  temperature,  and  the  polarization 
angle  of  the  light  in  the  fiber.  This  dependence  is  removed  by  summing  the  squares  of  each 
of  the  interferometric  signals  (with  the  dc  component  removed),  resulting  in: 

M^+v^+w'  =  (1.19) 

and  using  this  result  as  a  normalizing  division  factor,  the  ratio  of  Equation  1.18  to 
Equation  1.19  is  then  integrated  to  produce 

Unity  scale  factors  have  been  assumed  on  all  the  operators  (summers,  dividers,  etc.) 
in  presenting  the  demodulation  technique. 
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Following  on  from  the  background  of  interferometry  and  demodulation,  the  theory 
and  demonstration  of  multiplexing  interferometric  sensors  using  intensity  modulation  of 
the  source  is  introduced. 


S 


li.  MULTIPLEXING  OF  INTERFEROMETRIC  SENSORS 


The  proposed  multiplexing  scheme  uses  'Optical  Fiber  Subcarrier  Multiplexing' 
(OSCUM)  [Reference  1.  pages  544-5,  and  Reference  3,  pages  167-8]  although  this 
approach  has  not  previously  been  demonstrated  with  interferometric  sensors.  Reference  2 
(section  4.5,  pages  167-8)  details  the  use  of  this  multiplexing  technique  with  intensity 
sensors. 

A.  THEORY 

The  theory  will  be  presented  using  Figure  2. 1  and  applying  subscripts  to  the  equations 
which  relate  to  the  numbered  'position'  in  the  Block  Diagram. 


Figure  2.1.  Multiplexing/Demiiltiplexing  Block  Diagram:  Circled  numbers  relate 
to  subscripts  in  Equations  2.1  -  2.9. 
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The  laser  is  driven  with  an  AM  modulation  current  given  as. 


4,(/)  =  +  msin(to.r)]. 

and  the  resulting  light  intensity  from  the  laser  I,(t)  is  then: 

A^{t)  =  i4,[l +  OTsin(tu^/)]. 

The  electrical  signal  from  the  output  of  the  sensor,  at  point  2,  is  then  given  as 

=  [>4  +  5cos(^)][l  +  /nsin(n;,,/)] 

=  A  +  Bcos{^)  +  mA  sin(fi;„/)  +  mB  cos(^)sin(6>„/), 

where  the  notation  of  Reference  1 1  has  been  used.  Here  (j>  denotes  our  signal  of  interest 
and  includes  the  drift  and  coupler  terms,  given  as; 

^  =^O+^d„y,(O  +  ^c  =  ^,sin(6),t)+^j„^(/)  +  /(120°),  (2.4) 

where  /  =  0,  +1,  or  -1.  Using  the  trigonometric  identity: 

sin(a)cos(^  =  ^sin(a-/9)  +  -t-sin(a+^),  (2  5) 

the  signal  at  point  2  can  then  be  rewritten  as. 

v^{t)  =  A  +  Bcos{^)  +  mA  sin  i(o„{)  +  -^[sin  (o}^  r  -  ^)  +  sin  (o)^  t +  <!>)],  (2.6) 

where  the  gain  as  a  result  of  the  transimpedance  amplifier  and  fiber  losses  are  all 
encompassed  in  the  amplitudes  of  each  term.  An  example  of  a  computer  generated 
spectrum  of  the  signal  represented  by  Equation  2.6  is  shown  in  Figure  2.2.  The  peak  signal 
phase  modulation  ^5  =  ti  radians,  the  drift  phase  term  is  zero,  the  signal  frequency  is  1 
kHz  and  the  modulation  frequency  is  30  kHz,  and  the  intensity  modulation  depth  m  = 
0.75.  Figure  2.2  clearly  shows  the  baseband  and  modulated  portions  of  the  signal. 


(2.1) 

(2.2) 

(2.3) 
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Figure  2.2.  Frequency  spectrum  of  an  interferometric  output  phase  modulated  7C 
radians  peak,  at  1  kHz,  and  amplitude  modulated  at  30  kHz. 

The  signal  v^O),  given  by  Equation  2.6,  is  now  mixed  with  a  scaled  voltage 
proportional  to  the  laser  drive  current  modulating  signal.  This  mixing  signal  has  no  dc 
term  and  is  defined  as: 

Vo(/)  =  Csin(<y„/).  (2.7) 

Thus  the  signal  at  point  3  becomes: 

Vj{t)=AC  sin(<y„/)  +  5Csin(a),,/)cos(^) 

V  M .  O  ) 

+  mAC  sin' {o)j)  +  mBC sin^(<y„/)cos(^). 


which  can  be  rewritten  as: 
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(29) 


mAC  mBC  /a 

Vj(0  =  — -  +  cos(<^) 


BC 

+  AC  sia(a}j)  +  — [sin(<y,/  +  ^>)  +  sin(^y^/  -  ^)] 

+  ^d£.cos(2(U^/)  +  — ^— [cos(26>^/  -  ^)  +  cos(2<u^/  +  ^)], 


where  the  top  and  bottom  lines  of  Equation  2.9  represent  the  modulated  portion  of 
Equation  2.6  (or  Figure  2.2)  frequency  translated  to  baseband  and  twice  the  carrier 
frequency,  respectively.  The  second  line  of  Equation  2.9  represents  the  baseband  portion 
of  Equation  2.6  (or  Figure  2.2)  translated  to  the  carrier  frequency.  Using  our  previous 
example  (Figure  2.2)  and  mixing  with  a  signal  of  frequency  of  30  kHz  (C  =  1  from 
Equation  2.7),  the  spectrum  shown  in  Figure  2.3  is  obtained.  Figure  2.3  clearly  shows  the 
translation  caused  by  mixing  and  how  subsequent  recovery  of  our  signal  of  interest  by  low 
pass  filtering  can  be  achieved.  It  will  be  shown  later  that  this  mixing  or  heterodyning  is 
necessary  because  the  modulated  signals  will  overlap  in  the  baseband,  precluding  their 
recovery  by  simple  filtering. 

The  other  two  legs  (relatively  phase  shifted  by  120°  from  each  other)  are  similarly 
mi.xed  and  filtered  and  can  then  be  demodulated  using  the  passive  symmetric  demodulation 
scheme  previously  described. 
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Figure  2.3.  Spectrum  of  signal  of  Figure  2.2  mixed  with  a  signal,  phase  and 
frequency  locked  to  the  laser  modulating  signal  at  30  kHz. 


B.  DEMONSTRATION  WITH  OPTICAL  SYSTEM 

It  was  verified  that  the  ratio  of  components  in  the  baseband  was  the  same  as  the 
modulated  band  (neglecting  the  DC  and  carrier  components)  by  use  of  an  analog 
electronic  simulation  of  the  optical  components  and  signals.  This  indicated  that  the  static 
phasing  would  be  preferred  with  this  approach.  Appendix  C  describes  this 
demonstration  and  contains  all  the  significant  results  obtained.  The  essential  result  of  the 
demonstration  was  the  verification  of  the  computer  generated  simulation  of  Figure  2.2. 
Once  this  demonstration  was  successfully  completed,  a  two-laser,  two-sensor  optical 
demonstration  was  tested;  the  schematic  is  shown  in  Figure  2.4.  Two  lasers  are  amplitude 
modulated  at  two  different  frequencies  f„,  and  respectively.  The  two  sensors 
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inleirogated  by  these  lasers  are  then  coupled  together  before  photodetection, 
demultiplexing,  filtering  and  demodulation. 


Figure  2.4.  Scheinatic  for  Optical  Demonstration. 

A  computer  simulation  example  of  two  sensor  signals  and  two  modulation 
(subcarrier)  frequencies  is  shown  in  Figures  2.5  and  2.6.  They  follow  directly  from  the 
derivation  of  the  previous  section  but  highlight  the  potential  for  overlap  of  sideband 


signals  and  the  heterodyning  technique  used  to  recover  the  sigruil  of  interest  This  example 
represents  the  multiplexing  and  demultiplexing  of  the  interferometric  signals  from  one  of 
the  three  legs  shown  in  Figure  2.4.  Figure  2.5  shows  two  modulated  interferometric 
signals  with  the  following  characteristics;  subcarrier  modulation  frequency  =  40  kHz, 
sensor  signal  amplitude  =  l  .Stt  radians,  sensor  signal  frequency /,y  =  1  kHz,  bias  angle 
(pji  =  45°,  and  f„2  =  60  kHz,  (f>,2  =  radians,/,^  =  600  Hz,  =  0°.  The  amplitudes  and 
bias  angles  of  the  two  sensor  signals  were  specifically  chosen  to  be  different  to  highlight 
the  overlap  at  baseband. 


Figure  2.5.  Example  of  two  modulated  interferometric  signals. 

Figure  2.6  shows  the  signal  of  Figure  2.5  after  heterodyning  with  modulating  signal  at 
frequency  The  overlapped  baseband  signals  have  been  translated  to  f„,  (40  kHz) 
and  the  interferometric  signals  of  interest  can  be  recovered  by  lowpass  and  bandpass 
filtering.  The  recovered  interferometric  signals  from  each  of  the  three  legs  can  then  be 
demodulated.  The  filtering  must  provide  the  same  phase  response  to  each  of  the  three  legs 
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(ideally  constant  over  the  frequency  range  of  interest)  so  that  the  interferometric  signals 
can  be  accurately  demodulated. 


0  20  40  60  80  100 


Frequency  (kHz) 

Figure  2.6.  Signal  represented  in  Figure  2.5  heterodyned  with  modulating  signal 
at  frequency  (40  kHz). 
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1.  Experimental  Setup 

The  equipment  setup  for  the  optical  demonstration  is  shown  in  Figure  2.7. 
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The  driver  and  modulation  section  uses  the  Seastar  LD200  Ultra  Stable  Power 
Supply  to  drive  the  semiconductor  lasers.  This  driver  can  modulate  the  intensity  up  to  a 
frequency  of  about  1  MHz.  The  optical  power  output  starts  to  drop  off  at  about  800  kHz 
The  variable  voltage  divider  in  Figure  2.7  is  required  because  of  the  conflicting 
requirements  of  the  AD633  multiplier  and  the  external  modulation  input  of  the  LD200 
The  external  modulation  input  of  the  LD200  converts  a  voltage  input  to  a  laser  drive 
current  output.  The  ratio  of  the  driving  voltage  to  resulting  drive  current  is  10  mV  to 
1  mA  [Reference  12].  Thus  a  voltage  of  the  order  of  50  to  100  mVpp  is  required  to 
provide  the  required  drive  current  modulation  levels.  The  AD633  multiplier  chip,  on  the 
other  hand,  prefers  to  see  the  largest  signal  possible  (up  to  10  V);  so,  by  use  of  the  divider, 
a  phase  matched  signal  can  be  simultaneously  provided  to  the  LD200  (50-100  mVpp)  and 
theAD633  (10  Vpp). 

Figure  2.8  shows  the  detection,  transimpedance  amplifier,  and  amplification  for  one 
of  the  three  legs.  The  two  gain  stages  are  included  to  amplify  the  signal  before  mixing  so 
that  a  reasonable  level  signal  (comparable  to  the  laser  modulating  signal)  is  seen  by  the 
multiplier  in  the  mixing  stage.  The  amplifiers  remove  the  DC  component  of  the  signal  and 
each  provide  a  noninverted  gain  of  20.  The  20  kQ  resistor  in  the  amplifier  input  path  is 
required  to  provide  a  return  to  ground  for  the  very  small  input  current  [Reference  14, 
section  3.05,  page  95].  The  noise  up  to  this  stage  and  introduced  by  the  amplification  will 
also  be  amplified,  i.e.,  there  will  be  no  increase  in  signal  to  noise  ratio  (SNR).  This  was  not 
considered  an  issue  in  this  demonstration  of  the  multiplexing  technique  but  would  have  to 
be  addressed  in  an  actual  sensor  system  demonstration. 

Figure  2.9  shows  the  mixing  and  low  pass  filtering  circuit  where  THB292  is  an 
active  filter  integrated  into  a  single  package  [Reference  13,  page  1 18]  with  a  3  dB  cutoff 
at  22  kHz.  The  filter  has  a  frequency  response  as  shown  at  Appendix  D.  This  low  cost 
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integrated  package  ($12  85  each)  was  chosen  to  simplify  the  implementation.  In  a 
practical  system  the  cutofT  required  may  well  be  far  greater  than  22  kHz.  A  rule  of  thumb, 
developed  by  Dr.  David  A.  Brown  and  extensively  verified  by  experience,  gives  the 
interferometric  signal  bandwidth  as  a  function  of  signal  frequency  and  amplitude  as. 

A  =  5^./,.  (2.10) 

where  the  subscript  sb  denotes  sigiuficani  bandwidth.  This  estimate  is  considered 
conservative  but,  nevertheless,  provides  a  reasonable  estimate  of  the  maximum  frequency, 
signal  level,  or  combination  of  both  that  can  be  demodulated  by  the  system  if  we  designate 
the  cutoff  frequency  of  the  low  pass  filter  as  X*  • 
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Figure  2.8  Detection  and  Traiisimpedniice  Ampliner  Stage. 


For  example,  a  system  with  a  dynamic  range  of  100  dB  at  20  kHz  with  an  assumed 
minimum  detectable  signal  of  5  prad/^Hz  and  a  maximum  of  4%  THD  (a  typical  value 
using  Demodulator  Version  2.1  [Appendix  B])  would  tolerate  a  maximum  signal 
calculated  as  follows; 
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f 

20  log 

V 


/ffad/ VHz^ 
5  /ffad  /  >/Hz  ^ 


100  dB  re  Imrad  /  VHz  at  20  kHz 
0.5  rad. 


and  a  lowpass  filter  cutoff  above : 


(2.11) 


=  5x0.5x20x10’  =  60kHz.  (2.12) 

This  cutoff  frequency  then  allows  for  a  maximum  signal  at  100  Hz,  and  a  maximum  of  4% 
THD,  of 

^  =  A  =  =  100  radians.  (2  13) 

5/,  500 

or  a  dynamic  range  of 

=  20x10*^  =  146dBrel/irad/VHzatl00Hz,  (2.14) 

5x10-^ 

assuming  the  same  noise  floor. 

The  benchtop  interferometers  of  Figure  2.7  are  described  in  detail  at  Reference  9, 
pages  60-66  (see  also  Appendix  D).  The  purpose  of  the  benchtop  interferometer  is  to 
provide  a  signal  from  a  push-pull  optical  sensor  that  can  be  easily  adjusted  by  a  function 
generator  connected  to  the  benchtop  interferometer.  The  three  outputs  of  the  circuit  of 
Figure  2.9  are  connected  to  the  9207  resistor  of  the  analog  synunetric  demodulator, 
version  2.1,  described  in  detail  in  Appendix  B.  This  allows  the  transimpedance  amplifier  of 
the  demodulator  to  be  reconfigured  as  a  voltage  amplifier  and  saves  having  to  rewire  or 
reconstruct  the  demodulator  front  end,  depending  on  whether  or  not  the  demodulator  is 
being  used  with  the  demultiplexing  scheme  or  the  direct  demodulation  of  interferometric 
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sensor  signals.  This  stage  oF  amplification  is  still  used  to  ensure  that  the  three  signals  have 
the  same  ac  amplitude,  just  as  in  the  direct  demodulation  configuration  [Appendix  B]. 


V* 


LCTECreO  and 
AK'1_IFIED  OUTPUT 

LASER  MODXATBC 
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Figure  2.9.  Circuit  of  the  mixing  and  low  pass  filteriug  stage  of  Figure  2.7. 

2.  Results 

Prior  to  coupling  the  two  benchtop  interferometers  together,  each  was  tested 
separately.  A  modulation  signal  was  applied  to  the  appropriate  laser  and  then  the  sensor 
signal  was  demultiplexed  and  demodulated.  Figures  2. 10  and  2. 1 1  show  the  demultiplexed 
demodulated  outputs  for  inputs  to  both  channels  prior  to  coupling  together. 
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200  }.isec/div. 


Figure  2.10.  The  upper  signal  represents  the  input  to  benchtop  interreroineter  ^1 
(<|>,i  ==  n  radians  amplitude)  while  the  lower  signal  represents  the 
demultiplexed/demodulated  output  (f_,  -  80  kHz). 


500  nsec/div. 

Pigtire  2.1 1.  Upper  sigiinl  represents  the  input  to  beiichtop  interferometer  U1 
(({>,2  =  2.5  n  rndi.'ins  amplitude)  while  the  lower  signal  represents  the 
dcniultiplexed/deniodulated  output  (r„,2’=  120  kHz). 


After  connecting  the  two  benchtop  interferometers,  with  a  50/50  coupler,  the 
simulated  acoustic  signals  for  each  interferometer  were  again  demultiplexed  and 
demodulated.  The  outputs  for  two  different  inputs  to  each  of  the  benchtop  interferometers 
are  shown  in  Figures  2. 12  and  2. 13. 
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Figure  2.12.  Interferometers  1  and  2  coupled  together.  The  upper 
signal  represents  the  input  to  benchtop  interferometer  Hi  (4»,i  =  2.5  n 
radians  amplitude)  while  the  lower  signal  represents  the 
demultiplexed/demodulated  output  (f^,  =  80  kHz). 
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Figure  2. 14  shows  the  demultiplexed  and  demodulated  output  for  a  triangular  wave 
benchtop  interferometer  input  (2.5  V  at  -1350  Hz).  Even  though  the  frequency  spectrum 
of  triangle  wave  has  numerous  components  above  its  fundamental  (1350  Hz)  and  the  low 
pass  filter  cutoff  frequency  is  22  kHz,  the  demodulator  is  able  to  demodulate  the  band- 
limited  demultiplexed  system.  This  figure  was  chosen  to  highlight  the  robustness  but 
ultimate  limitations  of  the  technique. 

aii_  IV  nc  nuL  n  saoun  -577nv  fkti 
rM2  so««vnc 


Figure  2.14.  Interferometers  1  and  2  coupled  together.  The  upper  signal 
represents  the  input  to  benchtop  interferometer  ffl  while  the  lower 
signal  represents  the  demultiplexed/demodulated  output  (f^,  »  80  kHz). 
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Figures  2. IS  and  2.16  show  the  baseband  and  demultiplexed  signals,  respectively, 
for  a  500  Hz  benchtop  interferometer  signal.  The  ratios  of  the  components  of  the  signals 
appear  to  be  in  the  same  ratio  but  Table  2.1  shows  they  are  not.  The  reader  is  reminded 
that  this  was  not  the  case  when  the  multiplexing  scheme  was  tested  with  analog 
electronics.  This  difference  in  the  ratios  of  the  components  may  be  due  to  a  number  of 
causes  (or  a  combination  thereof)  These  include  the  phase  response  of  the  low-pass  filter 
not  being  constant  over  the  entire  passband,  the  optical  path  difference  (OPD)  in  the 
interferometer  legs  and  the  drift  angle  change  that  occurred  between  plots.  Independent  of 
the  reason  for  the  change  in  ratios  of  components,  the  change  is  the  same  in  all  three  legs 
of  the  interferometer  and  does  not  affect  the  ability  of  the  demodulator  to  correctly 
demodulate  the  signal.  This  would  tend  to  indicate  that  the  cause  is  due  to  the  OPD  in  the 
interferometer  or  the  change  in  drift  angle,  rather  than  the  phase  response  of  the  filters 
Further  investigation  of  the  cause  of  this  change  was  not  carried  out  since  it  was  always 
the  same  in  all  three  legs  and  essentially  amounted  to  a  drift  angle  change  between  the 
baseband  and  demultiplexed  signal. 

Although  Figures  2.15  and  2.16  do  not  show  the  noise  floor,  the  measured 
difference  between  noise  floors  in  the  baseband  spectrum  (-1 14  dB  re  1  V/VHz  at  5  kHz) 
and  the  demultiplexed  spectrum  (-98  dB  re  1  V/VHz  )  was  16  dB  at  5  kHz  (see  Appendix 
D,  Figures  D8  and  D9).  A  noise  floor  of  -98  dB  re  1  V/VHz  equates  to  minimum 
detectable  signal  of  13  ttrad/VHz.  This  is  a  large  difference  but  is  seen  to  be  due  to  the 
electronics  since  the  optical  crosstalk  between  sensors  was  less  than  60  dB.  The  crosstalk 
was  measured  by  mixing  with  f„i  and  measuring  the  size  of  the  known  spectral 
components  due  to  interferometer  #2.  The  electronics  for  this  demonstration  was 
assembled  on  a  bread  board  with  many  leads  to  and  from  the  demodulator.  This 
preliminary  measurement  of  the  noise  floor  is  considered  only  an  estimate  of  the  system 
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noise  floor  The  increase  in  noise  floor  and  hence  detectability  threshold  due  to  the 
nuiltiplexing  scheme  should  not  be  more  than  about  3-4  dB  as  a  result  of  the  modulation 
losses  and  demultiplexing  electronics. 
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Figure  2.15.  Baseband  speciriint  of  500  Hz  benchtop  interreroineter 
signal  (Note  60  IIz  noise)  prior  to  multiplexing,  =  0.6  radians. 
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Figure  2.16.  Spectrum  of  signal  in  Figure  2.15  after  demultiplexing. 
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Table  2.1  Frequency  component  ratios  for  Figures  2.15  and  2.16. 


C.  PRACTICAL  LIMITATIONS 

As  a  result  of  the  optical  demonstration  several  important  practical  limitations  of  the 
multiplexing  scheme  became  evident.  These  considerations  did  not  have  any  significant 
effect  on  the  demonstration  of  the  technique  but  would  be  important  considerations  in  a 
complete  system  demonstration.  The  most  important  of  these  is  semiconductor  laser  diode 
characteristics  and  the  subsequent  system  bandwidth  limitations. 
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1.  Semiconductor  Laser  Diode  Characteristics 
A  plot  of  the  optical  power  output  ver^s  input  drive  current  curve  for  an  ideal  laser 
diode  would  indicate  that  the  maximum  amplitude  modulation  and  corresponding  bias 
point  utilize  the  complete  linear  range  of  laser  operation  [Figure  2.17],  Figure  2.18  shows 
an  operating  curve  that  has  been  exaggerated  to  highlight  the  nonlinearities  and  power 
saturation  curvature. 


Figure  2.17.  Bias  point  and  amplitude  modulation  range  for  a 
laser  diode  (after  Reference  16). 


The  'kinks'  or  nonlinearities  in  the  operating  region  are  a  result  of  inhomogeneities  in 
the  active  region  and  power  switching  between  different  lateral  modes  of  the  laser.  The 
downward  curvature  at  the  upper  end  of  the  operating  curve,  or  power  saturation,  is 
attributed  to  active  layer  heating  [Reference  16,  page  105]. 
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Figure  2.18.  Example  of  kinks  and  power  saturation  in  laser  diode 
operating  curve  (after  Reference  16). 

Details  of  the  causes  and  methods  to  reduce  these  nonlinearities  is  well  documented 
(e  g.,  Reference  1,  chapter  5;  Reference  5,  chapter  3;  and  Reference  16,  chapter  4). 
Operating  curves  for  all  the  lasers  used  in  the  demonstration  are  included  (Appendix  D) 
but  are  far  too  crude  to  appreciate  the  significance  of  these  nonlinearities.  These 
nonlinearities  cause  higher  harmonics  of  the  modulation  frequency.  Modulation  of  sensor 
information  at  higher  harmonics  of  the  carrier  represents  lost  energy  that  cannot  be 
recovered  by  the  demultiplexing  and  demodulation  process. 

The  nonlinearities  can  be  better  quantified  by  measuring  the  frequency  spearum  of 
the  optical  output  of  an  intensity  modulated  laser  and  comparing  this  to  the  fi-equency 
spectrum  of  the  modulating  drive  current.  Figure  2.19  shows  the  frequency  spectrum  of 
the  HP3314A  function  generator  used  to  generate  the  intensity  modulation  of  the  lasers 
shown  in  Figures  2.20  and  2.21.  The  spectrum  of  the  HP3314A  is  within  specifications 
[Reference  17,  page  74], 
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The  two  lasers  shown  are  typical  of  semiconductor  laser  diodes  with  the  second 
harmonic  approximately  20-30  dB  down  from  the  fundamental  being  the  typical  range 
stated  in  several  of  the  references.  One  of  the  best  traditional  techtuques  to  reduce  this 
harmonic  distortion  to  30-40  dB  involved  the  circuit  technique  of  quasi-feedforward 
compensation  [Reference  18]  but  this  technique  has  been  superseded  with  the  advent  of 
special  purpose  highly  linear  devices.  The  linearity  of  such  devices  is  measured  in  terms  of 
the  third  harmonic  intercept  (TOI)  [see  Appendix  D,  Figure  D7].  These  lasers  are  more 
expensive  ($750-$1200  US)  but  may  be  required  depending  on  the  complete  system 
requirements.  It  should  also  be  noted  that  these  lasers  can  only  be  obtained  in  the  1300 
and  1550  nm  wavelengths. 

Also,  if  less  expensive  sources  are  used,  the  higher  order  harmonics  of  the  carrier 
impose  a  bandwidth  limitation  on  the  system.  Modulation  of  different  sources  has  to  be 
carried  out  in  one  of  two  ways;  either  at  frequencies  between  the  lowest  carrier  and  its 
second  harmonic,  or  at  carriers  separated  from  each  other  by  more  than  twice  the 
frequency  of  the  next  lower  carrier  frequency.  The  latter  method  can  be  used  since  the 
third  and  higher  harmonics  are  generally  negligible.  Both  these  methods  impose  a  stricter 
limit  on  the  number  of  possible  sensors  than  with  the  use  of  highly  linear  sources  where 
higher  order  carrier  components  are  negligible. 
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Figure  2.19.  Frequency  Spectrum  orHP3314A  Function  Generator  for 
desired  output  of  10  Vpp  at  100  kllz  used  to  drive  lasers  of 
Figures  2.20  and  2.21. 
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Figure  2.20.  Frequency  SpecCrum  orififaclii  IfL83l2G  #SP27f2  Inser 
driven  by  Seastnr  LD200  power  supply  fed  by  signal  of 
Figure  2.19. 
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Figure  2.21.  Frequency  Spectrum  of  Sharp  LT015MDO  ttSP2676  laser 
driven  by  Seastar  L0200  power  supply  fed  by  signal  of 
Figure  2.19. 
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2.  Aging  and  Coherence  Length  Considerations 
The  coherence  length  of  the  laser  must  be  greater  than  the  optical  path  difference  in 
the  interferometer  legs.  Using  the  typical  data  shown  at  Reference  19,  page  6S.  a 
coherence  length  of  approximately  7  mm  could  be  inferred.  This  was  calculated  using  a 
full  line  width  dA  of  0.1  nm  (estimated  from  Reference  19)  and  an  operating  wavelength  of 
830  nm  from; 


This  Lq  of  7  mm  was  presumed  to  be  close  to  the  optical  path  difference  of  the  benchtop 
interferometer.  The  causes  and  effects  of  aging  of  semiconductor  laser  diodes  are  well 
documented  [Reference  1,  Volume  I,  section  5. 3. 3. 4,  Reference  5,  page  56,  and  Reference 
16,  section  4-5].  One  cause  is  related  to  the  drive  current  level.  The  laser  had  been 
operated  at  high  drive  levels  and  this  may  have  possibly  reduced  the  coherence  length  of 
the  laser.  The  laser  was  replaced  with  a  Hitachi  HL8312G  which  was  presumed  to  have  a 
longer  coherence  length  (i.e.,  greater  than  the  OPD  in  the  interferometer),  if  only  because 
it  was  new.  The  fringe  visibility  returned.  This  problem  highlighted  the  importance  of 
coherence  length  and  optical  path  difference  considerations  and,  especially,  their 
relationship  as  a  function  of  aging  due  to  high  drive  current  levels,  high  temperature 
operation,  or  any  other  aging  mechanism  which  will  change  the  operating  characteristics 
of  the  laser  diode. 

3.  Operating  Wavelength  Considerations. 

The  choice  of  830  nm  as  the  operating  wavelength  for  the  optical  demonstration  of 
this  thesis  was  motivated  by  the  inexpensive  CD  type  semiconductor  laser  sources 
available  at  this  wavelength.  It  is,  however,  very  important  to  consider  the  overall  cost  of 
the  various  system  components  and  the  direction  in  which  the  optical  fiber 
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communications  industry  and  technology  are  moving.  Fiber  for  1300  nm  sources  is  now 
cheaper  than  that  for  830  nm.  The  highly  linear  source  described  earlier  are  only  available 
in  the  1300  and  1550  nm  ranges.  Optical  amplifiers,  which  may  be  a  consideration  in 
multiplexing  (considered  later),  are  presently  only  available  in  the  1550  nm  range. 

All  these  cost  and  component  availability  issues  need  to  be  considered  when 
choosing  the  operating  wavelength  for  an  all  optical  sensor  array.  It  may  well  be  that  the 
developments  of  the  communications  industiy  will  determine  the  most  economical 
operating  wavelength.  In  either  case  this  multiplexing  technique  provides  a  less  expensive 
source  option.  This  technique  indicates  the  possibility  for  use  with  inexpensive  CD  type 
sources  at  830  nm  ($150-350  US),  or  relatively  inexpensive,  highly  linear  sources  similar 
to  those  for  cable  television  applications  at  1300  or  1550  nm  ($750-1500  US).  These 
sources  are  both  inexpensive  compared  to  the  highly  coherent  sources  required  by  other 
presently  used  PGC  FDM  techniques  involving  wavelength  modulation  of  the  source  [e  g.. 
Reference  2,  pages  170-172],  in  which  sources  typically  cost  $10,000-15,000  US.  This 
cost  saving  is  significant  enough  to  warrant  further  research  into  applying  the  technique  to 
larger  arrays  at  RF  frequencies.  With  the  present  cost  of  the  sources,  fiber,  demodulator 
and  demultiplexing,  the  most  expensive  component  in  the  system  is  the  3x3-coupler  in 
each  sensor  which  presently  cost  $400-750  US  each,  while  2x2-couplers  cost  on  the  order 
of  $100  US. 


38 


in.  DESIGN  OF  AN  M  X  N  ARRAY 


A.  MULTIPLEXING  ARCHITECTURE  FOR  AN  M  x  N  ARRAY 

The  multiplexing  technique  demonstrated  optically  and  described  in  Section  II  B, 
utilizing  two  lasers  and  two  sensors,  can  be  extended  to  an  A/  x  ^  sensor  element  array. 
Figure  3.1  shows  a  possible  implementation  topology  for  a  3x4  element  array.  Each  of 
the  bold  vertical  lines  connecting  a  stave  of  sensors  (or  vertical  column  in  Figure  3.1) 
actually  represents  three  separate  fibers  and  has  only  been  drawn  as  one  for  clarity  in  the 
figure.  The  three  interferometric  signals  obtained  for  each  stave  contain  the  sensor 
information  from  all  the  sensors  in  that  stave.  To  reconstruct  the  individual  sensor 
information  the  three  interferometric  lines  must  be  demultiplexed  once  and  low-pass  or 
band-pass  filtered  as  appropriate  before  demodulation.  Alternatively  each  sensor's 
interferometric  signals  could  be  demultiplexed  and  low-pass  filtered  prior  to 
demodulation.  For  either  approach  to  demultiplexing,  a  demodulator  is  required  for  each 
sensor  in  the  array  if  ail  are  to  be  required  to  be  interrogated  simultaneously. 

This  is  not  considered  to  detract  from  the  technique  since  the  cost  of  parts  of  each 
demodulator  [Appendix  C]  is  approximately  $150  US  [Reference  21]  and  can  be  made 
cheaper  with  dedicated  integrated  circuits.  This  is  a  significant  cost  reduction  over  an 
earlier  version  of  the  demodulator  [Reference  9,  page  204]  which  quoted  the  cost  as  $270 
US.  This  cost  could  possibly  be  reduced  further  by  several  means,  including  large  quantity 
manufacture,  special  purpose  digital  signal  processing  (DSP)  cards,  or  special  purpose 
DSP  chips.  The  approach  of  a  specific  DSP  card  follows  on  from  the  Naval  Postgraduate 
School  research  conducted  by  Lt.  Brian  R.  McGinnis,  USN,  and  is  presently  being 
undertaken  as  a  thesis  project  by  LCDR  David  W.  Brenner,  Canadian  Armed  Forces. 


Figure  3.1  Multiplexing  Architecture  for  a  3  x  4  Sensor  Array. 

The  possibility  of  reducing  the  number  of  demodulators  by  a  system  of  time  division 
multiplexing  (TDM)  of  sensor  information  is  also  an  option.  The  added  complexity  and 
cost  of  incorporating  TDM  may  be  greater  than  that  of  the  extra  demodulators. 
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Similarly,  the  relatively  low  cost  of  optical  fiber  for  short  length  systems,  i  e  ,  those 
where  sensor,  source,  and  receiver  separation  is  of  the  order  of  tens  of  meters,  has  not 
made  the  reduction  in  the  number  of  return  lines  a  critical  design  parameter.  There  may  be, 
however,  an  operational  requirement  to  reduce  the  number  of  return  lines.  A  deployable 
array  or  remotely  operated  vehicle  (ROV)  transmitting  data  to  a  remote  location  are 
possible  examples  where  such  a  design  consideration  may  be  essential  Methods  to  achieve 
significant  reduction  in  return  lines  are  considered  later. 

An  analysis  of  the  topology  of  the  architecture  proposed  in  Figure  3.1  can  be 
conducted  to  demonstrate  the  savings  and  costs  associated  with  the  multiplexing 
technique.  The  ratio  of  sensors  to  channels  (number  of  optical  fibers)  represents  a  measure 
of  the  saving  due  to  multiplexing  that  can  easily  be  quantified.  Using  the  architecture  of 
Figure  3. 1,  the  number  of  lasers  is  denoted  by  M  and  the  number  of  sensors  per  laser  by  N. 
The  number  of  fiber  lines  F  is  then  given  as  the  number  of  lasers  plus  three  times  the 
number  of  staves  or  return  lines; 


F  =  M  +  3A^. 


(3  1) 


Thus  for  the  3  X  4  array  of  Figure  3.1  we  find  the  number  of  sensors  S  and  number  of 
fiber  lines  as; 


and 


s  = 

MxN  =  12, 

(3.2) 

F  = 

M-i-3N  =  15, 

(3.3) 

MxN 

1  y  .  a  r 

=  0.8  sensors  /  channel . 

(3.4) 

F  M  +  3N 
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For  a  10  X  10  array; 


s  = 

MxN  =  100, 

(3.5) 

F  = 

M  +  3N  =  40, 

(3.6) 

100 

A/\ 

=  2.5  sensors  /  chaimel . 

(3.7) 

Thus  as  the  number  of  sensors  increases  the  sensor-to-channel  ratio  increases  and  a 
multiplexing  gain  is  realized. 


Table  3.1  Number  of  fiber  lines,  F,  for  M  x  A'  array 


1  Sensors/laser  Number  of  lasers  M  \ 

N 

2 

3 

4 

5 

6 

7 

8 

9 

10 

2 

8 

9 

10 

11 

12 

13 

14 

15 

16 

3 

11 

12 

13 

14 

15 

16 

17 

18 

19 

4 

14 

15 

16 

17 

18 

19 

20 

21 

22 

5 

17 

18 

19 

20 

21 

22 

23 

24 

25 

6 

20 

21 

22 

23 

24 

25 

26 

27 

28 

7 

23 

24 

25 

26 

27 

28 

29 

30 

31 

8 

26 

27 

28 

29 

30 

31 

32 

33 

34 

9 

29 

30 

31 

32 

33 

34 

35 

36 

37 

10 

32 

33 

34 

35 

36 

37 

38 

39 

40 

It  should  also  be  noted  that  this  system  would  work  with  the  NRL  quadrature 
demodulation  approach  [Reference  6]  which  uses  two  of  the  three  lines. 
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1.  Power  Budget  Requirements. 

The  limitation  on  the  number  of  sensors  possible  is  a  function  of  the  laser  power 
available,  the  number  of  couplers  required,  and  the  bandwidth  or  frequency  separation 
between  modulation  frequencies.  The  most  efficient  power  split  ratio  for  each  coupler  can 
be  determined  from  its  location  in  the  array.  The  coupler  ratio  required  for  each  coupler  in 
a  chain  (horizontal  row  in  Figure  3.1)  from  any  laser  can  be  calculated  using  Figure  3.2  as 
follows  (neglecting  intrinsic  losses); 


A»)  = 

N 


(3.8) 


where  P(n)  is  the  power  transmitted  by  coupler  n  to  the  next  coupler,  the  output 
power  of  the  laser,  and  P^  the  (constant)  power  transmitted  to  each  sensor  The  coupler 
ratio  C„  is  then  defined  as  the  ratio  of  the  transmitted  power.  P(n),  to  the  input  power, 
P(n-l),  and  is  given  by. 


C. 


PW  ^  t  P. 
P(n-\)  P(n-\) 


— — ; -  for  //  =  0,1,  ,N-\. 

+ 1  -  // 


(3.9) 


Normalizing  these  calculations  for  a  given  laser  output  power  of  P„  .  Equation  3.8 
becomes; 


F(//) 


P{n) 

Po 


(3-10) 


where  P’(n)  in  Equation  3.10  represents  the  fraction  of  the  total  laser  power  available  at 
coupler  n. 
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Figure  3.2  Power  Transmission  Curve  for  Equal  Power  to  Each  Coupler. 

The  calculations  for  the  coupler  power  split  ratios  in  of  the  return  lines  of  each  stave 
(vertical  column  of  sensors  in  Figure  3.1)  can  be  calculated  by  using  the  formulas  from  the 
previous  calculations  for  the  coupler  ratios  in  a  laser  chain  but  for  M  sensors.  These 
coupler  ratios  must  then  be  implemented  in  reverse  order  down  a  stave.  Thus  the  first  or 
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top  couplers  in  a  stave  of  sensors  will  have  a  coupling  ratio  of  0.5  while  the  last  or  bottom 
couplers  will  have  a  coupler  ratio  of  1/M . 

In  a  particular  laser  chain,  power  is  not  wasted  in  the  sense  that  light  is  either 
coupled  to  a  sensor  or  transmitted  to  sensors  ftirther  along  the  chain.  The  only  losses  in 
the  chain  are  splice,  coupler,  and  fiber  losses.  In  contrast  the  couplers  in  a  stave  each  have 
one  leg  of  wasted  power.  The  most  efficient  use  of  optical  power  therefore  will  maximize 
N  and  minimize  M.  The  compromise  is  that  the  sensor-to-channel  ratio  is  reduced  for 
larger  N  and  smaller  M.  A  compromise  between  the  two  design  parameters  of  maximum 
optical  power  per  sensor  and  maximum  sensor*to-channel  ratio  must  be  made  based  on  the 
application  requirements. 

A  MATLAB'''^  routine  [Appendix  D]  was  written  to  calculate  the  maximum  output 
power  from  each  sensor  for  any  given  array  size,  M  x  N.  This  routine  takes  into  account 
the  system  losses  based  on  the  following  estimates  of  the  individual  losses.  These 
estimates  are  based  on  the  typical  losses  for  those  components  used  in  the  2  x  1  optical 
demonstration.  Loss  estimates  are  assumed  as  follows: 

0.5  dB 
0.2  dB 
3.0  dB/km 
80  m  [Reference  22] 

3  mW  or  4.8  dBm 
-1 10  dBV  re  1  prad/VHz 
30  m  (e  g.,  for  50  Hz  array) 


•  loss/coupler  (single  mode) 

•  loss/splice  (single  mode  fiber) 

•  fiber  loss  @  X  =  830  nm 

•  fiber  length/sensor 

•  typical  laser  output  power(at  Im  pigtail) 

•  typical  noise  floor  at  receiver 

•  fiber  length  between  sensors 

Thus  the  loss  per  sensor  can  be  calculated  as: 


2 x0.5(coupler)  +  4 x0.2(splice)  +  3 x0,08km(sensor fiber)  =  2.04 dB/ sensor,  (3.1 1) 
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and  the  loss  per  coupler  including  the  fiber  loss  between  couplers  as: 

0.5  +  3  X 0.03 km( fiber  loss)  =  0.59  dB.  (3  12) 

The  loss  per  coupler  (0.59  dB)  and  loss  per  sensor  (2.04  dB)  have  been  included  in  the 
MATLAB™  routine. 

Examples  are  shown  below  for  different  values  of  M  and  N  for  both  the  fraction  of 
laser  power  available  and  an  estimated  power  of  3  mW,  based  on  data  for  the  HL8312G 
Hitachi  source  [Reference  23].  By  comparing  the  calculations  for  a  4  x  25  element  array 
and  a  10  x  10  element  array,  the  compromise  between  optical  power  and  number  of  return 
lines  for  these  two  implementations  of  a  100  element  array  can  be  seen  Figures  3  3  and 
3.4  compare  the  fractional  power  transmitted  to  a  particular  coupler  n  in  a  laser  chain  for 
the  two  cases.  Figures  3.5  and  3.6  compare  the  coupler  ratios  for  couplers  in  a  particular 
laser  chain  for  the  two  arrays.  The  reason  that  the  last  coupler  is  not  exactly  one  half  in 
each  case  is  due  to  the  estimated  losses  included  in  the  calculations.  Figures  3.7  and  3.8 
show  the  power  from  each  sensor  in  a  laser  chain  for  the  two  cases.  There  is 
approximately  2.7  dBm  more  power  delivered  to  each  sensor  in  the  case  for  A^=  10. 
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Fractional  Power  Input  Pin/Rotal 
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Figure  3.7.  Estimated  power  coupled  to  each  sensor  in  chain  for  Hitachi 
HL8312G  laser  (3  mW  at  1  m  pigtail),  including  estimated  losses,  N  =  25. 
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Figure  3.8.  Estimated  power  coupled  to  each  sensor  in  chain  for  Hitachi 
HL8312G  laser  (3  mW  at  1  m  pigtail),  including  loss  estimation,  N  =  10. 


Figures  3.9  and  3.10  show  the  coupler  ratios  for  each  coupler  down  a  stave  (or 
column  in  Figure  3.1).  Again  the  first  coupler  is  not  exactly  equal  to  one  half  as  a  result  of 
the  losses.  The  most  important  results  of  this  analysis  are  shown  in  Figures  3  11  and  3  12 
where  the  estimated  optical  power  for  each  sensor  down  a  stave  are  compared  for  the  two 
cases.  The  power  for  the  last  sensor  in  the  stave  represents  the  power  coupled  from  each 
leg  of  each  sensor  in  the  array  that  is  received  at  the  photodetector  It  can  be  seen  that  for 
the  10  X  10  array  the  power  lost  as  a  result  of  the  lost  (unused)  power  from  each  coupler 
in  every  stave  is  so  significant  that  there  wilt  insufficient  power  at  the  photodetector.  The 
estimated  system  noise  floor  is  expected  to  be  approximately  -90  to  -1 10  dBm  The  4  x 
25  array,  on  the  other  hand,  is  able  to  deliver  approximately  -58  dBm,  which  is  on  the 
order  of  50  dB  above  that  for  the  10x10  array. 
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Figure  3.9.  Coupler  ratio  for  couplers  in  each  stave  (column  in 
Figure  3.1)  for  4  x  25  element  array. 
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Figure  3.10.  Coupler  ratio  for  couplers  in  each  stave  (column  in 
Figure  3.1)  for  10  x  10  element  array. 


An  optical  sensor  array  topology  is  not  only  dilTerent  from  present  piezoelectric  with 
twisted  pair,  or  coaxial  lines,  but  may  vary  depending  on  the  specific  application.  A  hull 
mounted  array,  for  example,  may  best  be  implemented  with  no  multiplexing  in  either  one 
laser  per  sensor,  or  a  shared  laser  topology  as  shown  in  Figure  4.1.  The  savings  in  size  and 
weight  realized  by  an  optical  system  over  those  presently  available  may  preclude  the 
requirement  for  multiplexing 


Figure  3.13.  Array  topologies  for  a  hull  mounted  optical  array: 
one  sensor  per  laser  (upper)  and  shared  laser  chain  (lower). 
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2.  Bandwidth  Limitations 

Consider  a  system  of  M  lasers  with  the  lowest  modulation  frequency  j  at  20 
MHz.  Let  the  upper  design  limit  of  20  kHz  be  the  maximum  acoustic  frequency  with  a 
dynamic  range  of  140  dB  at  100  Hz  Using  Equation  2.10,  and  assuming  similar  noise 
floor  and  maximum  acceptable  distortion  levels,  this  gives  the  significant  bandwidth  of  the 
interferometric  signal  as; 

A  =  =  5x(0.5x10^x10““'“)x20  =  500  kHz  (3.13) 

With  the  fact  that  we  have  both  upper  and  lower  sidebands,  a  1  MHz  frequency  separation 
is  required  between  laser  modulation  frequencies.  For  at  20  MHz  this  represents  an 
upper  limit  of  20  sensors  that  can  be  accommodated  before  the  second  harmonic  of  the 
lowest  carrier  . 

Investigation  of  an  M  x  array  reveals  that  the  optical  power  restrictions  imposed 
by  coupling  sensors  in  a  stave  and  the  subsequent  lost  or  wasted  power  at  every  coupler  is 
a  far  more  stringent  design  consideration  than  bandwidth  limitations 

B.  REDUCTION  IN  THE  NUMBER  OF  RETURN  LINES. 

An  application  such  as  a  deployable  array  or  remotely  operated  vehicle  (ROV)  may 
require  significant  reduction  in  the  number  of  return  lines.  The  number  of  return  lines  may 
be  required  to  be  reduced  beyond  that  shown  in  Figure  3. 1  or  even  down  to  a  single  return 
line.  The  cost  of  additional  optical  fiber  may  be  greater  than  that  of  the  extra  couplers  and 
processing  required  in  a  system  with  reduced  return  lines.  This  is  the  case  in  systems 
where  the  distance  between  sensors  and/or  receivers  and/or  sources  become  too  large.  At 
approximately  50  mik  iv  the  gallon  (container),  fiber  is  far  cheaper  and  more  co'np.'ct 
than  copper. 
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A  simple  method  for  reducing  the  number  of  return  lines  utilizing  3x3-coupler 
terminated  sensors  and  symmetric  demodulation  is  a  combination  of  OSCUM  detailed 
earlier  and  time  division  multiplexing  (TDM).  TDM  has  been  successfully  demonstrated 
with  intensity  based  sensors  [e  g..  Reference  2,  pages  165-6,  and  Reference  4,  page  95] 
and  interferometric  sensors  [e  g..  Reference  2,  page  172-176,  Reference  3,  pages  82-84, 
and  Reference  4,  pages  96-100]. 

The  discussion  of  Reference  3,  page  83,  in  which  TDM  with  Mach-Zehnder  sensors 
terminated  with  2x2-couplers  is  presented,  can  be  adapted  to  the  case  for  3x3-coupler 
terminated  sensors  and  symmetric  demodulation.  Such  a  scheme  is  shown  in  Figure  3.14. 
The  demultiplexing  of  the  pulses  would  have  to  be  carried  out  prior  to  demodulation 
Each  input  pulse  generates  one  pulse  per  sensor  on  each  of  the  three  output  lines.  These 
output  pulses  must  be  appropriately  gated  to  relate  them  to  their  respective  sensor  In 
order  to  correctly  determine  which  pulse  is  associated  with  which  sensor  and  to  ensure 
that  there  is  no  overlap  of  pulses,  the  pulse  duration  must  be  less  than  the  propagation 
delay  between  any  two  sensors,  and  the  time  between  pulses  must  be  greater  than  the 
propagation  delay  difference  between  the  shortest  and  longest  paths. 

The  power  output  per  sensor  can  be  maximized  in  a  fashion  similar  to  the  previous 
section  by  choosing  the  coupling  ratios  of  the  couplers  appropriately.  There  is,  however, 
an  additional  power  loss  in  a  TDM  system  due  to  the  duty  cycle  of  the  source. 
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Figure  3.14.  TDM  willi  3x3-cou|}ler  terminated  sensors  for  symmetric 
demodulation  (adapted  from  Reference  3). 


1'he  number  of  return  lines  of  the  complete  array  could  be  reduced  to  a  single  line  by 
1  DM  of  the  individual  sensor  legs.  This  introduces  further  losses  and  demultiplexing 
complexity.  A  combination  of  OSCUM  and  TDM  could  be  employed  to  improve  the 
optical  power  per  sensor  (equivalently  the  SNR  of  each  sensor).  A  scheme  incorporating  a 
single  return  line  and  a  combination  of  OSCUM  and  TDM  is  proposed  in  Figure  3.15.  The 
loop  diameters  indicate  relative  differences  in  path  length  delays. 


n _ n _ TL. 


V/u vt'lei  iijU  I  (Jenu  - 
11  11  1  dti  iu  -; 

I  l€i'iO'  Jt  iki  Uv  u  I 


rigiire  3.15.  Single  rcluni  arcliitccliire  iiicorporaling  OFSUM  and  I'DM. 

It  sliould  also  be  noted  that  a  system  with  a  reduced  number  of  return  lines  is  not  as 
robust  in  the  sense  that  the  loss  of  a  single  fiber  or  coupler  leads  to  a  much  larger  or  total 
loss  of  sensor  information  while  this  loss  is  reduced  for  the  case  of  more  return  lines. 
Again  the  application  will  dictate  the  costs  and  benefits  involved.  References  3,  4.  and  5 
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provide  an  overview  of  the  techniques  presently  used  with  a  variety  of  other  types  of 
sensors  and  demodulation  techniques.  Further  research  is  required  to  apply  these 
techniques  to  3x3-coupler  terminated  Mach-Zehnder  sensors  with  symmetric 
demodulation. 
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rv.  CONCLUSIONS 


The  technique  of  'optical  fiber  subcarrier  frequency  multiplexing'  has  been  successfully 
demonstrated  with  a  2  x  1  element  array.  This  is  the  first  time  that  this  multiplexing 
technique  has  been  known  to  be  demonstrated  with  interferometric  sensors. 

The  architecture  for  a  particular  array  is  very  application  specific  with  the  requirements 
of  the  particular  application  dictating  the  multiplexing  technique  to  be  utilized.  This  occurs 
as  a  result  of  the  large  power  budget  trade-offs  incurred  when  optical  signals  are 
recombined  through  couplers.  The  development  and  introduction  of  low  cost  optical 
amplifiers  may  provide  a  solution  to  this  problem. 

The  following  list  details  some  of  the  follow-on  work  recommended  to  further  the  goal 
of  all  optical  sensor  systems  for  a  wide  variety  of  applications. 

•  Demonstration  of  the  technique  at  RF  frequencies  with  a  larger  number  of 
sensors  (or  simulated  sensors),  possibly  with  highly  linear  sot  a  as  similar  to 
those  described  in  Chapter  II,  section  B.2. 

•  Investigation  of  a  method  to  digitally  implement  the  multiplexing  scheme  to 
further  reduce  costs  and  hardware  complexity. 

•  Investigation  and  demonstration  of  techniques  to  reduce  the  number  of 
return  lines  for  specific  applications. 

•  A  less  expensive  3x3-coupler,  or  alternative  to  the  3x3-coupler  that 
presently  terminates  the  sensors  developed  at  NFS. 
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APPENDIX  A.  SPECTRUM  OF  AN  INTERFEROMETRIC  SIGNAL 


The  interferometric  signal  output  from  a  3x3  terminated  Mach-Zehnder  interferometer 
with  input  signal  <f)(t)  (Equation  Al),  expanded  into  a  Bessel  weighted  series,  is  given  by 
equation  A2; 


I(t)  =  A  +  Bcos(^t)) 
^t)  =  ^,sin<u,t. 


(Al) 


Joi^s)  +  2Y,JM)cosi2lco^t) 


/=! 


21-/  2/»l  (^,)sin  [(2/  +  l)<y,/] 


»=o 
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where  ^jri/i(0  is  the  slowly  varying  (compared  to  the  signal  of  interest)  or  drift  phase  term 
which  includes  the  phase  shift  due  to  environmental  effects  (hydrostatic  pressure, 
temperature,  etc.)  as  well  as  the  physical  difference  in  the  length  of  the  legs  and  the  phase 
shift  due  to  the  coupling  region  of  the  interferometer. 

Thus,  for  a  given  (assumed  to  be  a  constant  over  the  time  scale  of  the  signal 

of  interest  27t/<y^ )  the  ratio  of  the  spectral  components  of  odd  or  even  harmonics  of  the 
signal  of  interest  will  be  determined  by  the  ratio  of  the  Bessel  functions.  This  property  can 
be  used  to  determine  the  operating  point,  i.e.,  the  magnitude  of  by  measuring  these 
ratios  of  components  and  relating  them  to  the  equivalent  Bessel  function  ratios  (i.e.,  J2/J4, 
J3/J5,  etc.)  [Reference  9,  page  28].  This  is  especially  useful  for  determining  the  operating 
point  in  the  'sub-fnnge'  range.  Figure  A I  shows  the  lower  order  Bessel  functions  and 
Table  Al,  the  component  ratios  for  =  n  radians  amplitude,  which  can  be  used  to  verify 
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the  examples  of  *  Ji  radians  amplitude  and  drift  angles  ^dnfi  ~  120®,  and  240®, 

shown  in  Figure  A2. 


APPENDIX  B.  SYMMETRIC  ANALOG  DEMODULATOR  VERSION  2.1 


Reference  9  details  the  first  implementation  of  the  symmetric  demodulation  scheme 
using  analog  electronics.  The  analog  implementation  was  redesigned  by  Dr.  David  L 
Gardner  with  vastly  improves  the  small-signal  performance  over  the  original 
implementation.  Gardner's  implementation  was  not  previously  documented  and  was  used 
for  the  multiplexing  demonstrations  of  this  thesis. 

The  circuit  schematic  is  shown  at  Figure  Bl.  The  algorithm  is  well  documented 
[References  3,  4,  5,  and  7]  with  version  2.1  described  here  being  configured  for  lower 
noise  floor  and  ease  of  component  manipulation.  The  circuit  description  is  divided  into 
several  distinct  stages  adapted  from  Dr.  Gardner's  notes.  Performance  results  of  version 
2.1  are  included. 

A.  INPUT  STAGE 
1.  Photodiodes. 

Three  input  choices  are  available  as  follows. 

1.  Three  1300  nm  photodiodes  are  the  lower  row  of  bulkhead  SMA  connectors 
Presently  fitted  are  model  UDT080  devices  (with  an  active  area  of  0.08  mm^).  These 
devices  are  selected  by  opening  the  box  and  moving  each  of  the  three  blue  jumpers 
(adjacent  to  the  edge  connector)  to  the  left-most  position  as  one  looks  down  on  the  board 
facing  the  edge  connector.  These  small  and,  hence,  fast  diodes  may  present  some  difficulty 
in  alignment.  Both  the  fiber  and  diodes  can  be  epoxied  into  position  for  field  tests. 
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2.  Three  830  nin  photodiodes  are  the  upper  row  of  bulkhead  SMA  connectors 
Presently  fitted  are  model  HRS  102  devices  (with  an  active  area  diameter  of  300  pm) 
These  detectors  are  selected  by  opening  the  box  and  moving  each  of  the  three  blue 
jumpers  (adjacent  to  the  edge  connector)  to  the  right-most  position  as  one  looks  down  on 
the  board  facing  the  edge  connector.  These  also  should  be  fixed  into  position  with  epoxy 
for  field  tests. 

3.  It  is  often  convenient  to  use  interferometric  sensor  simulators  to  test 
demodulators  and  conduct  various  other  trials.  Since  the  output  fi-om  these  simulators  is  a 
voltage,  an  input  resistor  has  been  installed  to  allow  the  demodulator  to  be  used  with  the 
simulators.  This  will  cause  the  input  amplifier  to  be  an  inverting  type  rather  than  the 
transimpedance  configuration  required  for  the  photodiodes  Three  9207  Q  resistors  are 
precision  type  and  can  be  used  in  conjunction  with  the  fixed  value  feedback  resistor  to  test 
the  operation  of  the  demodulator  input.  The  box  must  be  opened  to  gain  access  to  these 
resistors.  They  are  gray  and  installed  vertically  along  the  edge  of  the  board  adjacent  to  the 
connector.  The  long  lead  of  the  resistor  is  available  for  connection  but,  first,  the  blue 
jumpers  must  be  removed. 

2.  Input  Amplifiers. 

Two  options  are  available  for  gain  control  of  the  input  amplifiers: 

1 .  Fixed-value  resistors  -  In  the  course  of  testing  interferometric  sensors,  it  is 
useful  to  determine  the  light  power  in  each  of  the  sensor  fibers.  The  typical  responsivity  of 
the  photodiodes  is  0.5  AAV.  Since  the  current-to-voltage  ratio  or  conversion  factor  must 
be  known  to  make  such  a  measurement,  20  kn  resistors  have  been  installed.  These 
resistors  are  selected  by  opening  the  box  and  moving  the  black  jumpers,  located  along  the 
left  edge  of  the  board  as  viewed  facing  the  edge  connector,  to  the  leftmost  position. 
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2.  Potentiometer-controlled  gain  -  The  output  power  from  each  leg  of  a  coupler 
is  rarely  uniform.  Variations  in  the  modulation  depth,  splice  losses,  bending  losses  and 
modulation  losses  all  reduce  the  light  emitted  from  each  fiber.  Since  the  demodulator 
algorithm  assumes  equal  power  from  each  leg,  potentiometers  have  been  installed  to 
provide  adjustment  of  the  amplitude  of  the  AC  portion  of  each  input  signal.  Adjustments 
of  the  DC  portion  may  also  be  accomplished.  The  potentiometers  are  blue  in  color  and  are 
located  just  in  from  the  left-hand  edge  of  the  board  as  viewed  facing  the  edge  connector. 
The  output  of  each  amplifier  may  be  monitored  by  connecting  an  oscilloscope  lead  to  the 
test  point  loops  near  the  yellow  jumpers  and  aligned  with  the  transverse  axis  of  the  board. 
These  test  points  are  connected  to  pin  5  of  the  non-inverting  amplifiers  (ICs  U4-U6  of 
Figure  Bl). 

B.  INTERMEDIATE  AMPLIFIERS  AND  OUTPUT  BUFFERS 

1.  Intermediate  Amplifiers 

The  intermediate  amplifiers  have  been  included  to  permit  adjustment  of  the  signal 
level  to  the  remainder  of  the  circuit.  The  gains  of  these  non-inverting  amplifiers  are 
selected  via  the  yellow  jumpers  towards  the  left  edge  of  the  board,  as  viewed  facing  the 
edge  connector.  The  gains  permitted  are  20,  30,  and  40  dB,  using  pairs  of  jumper  leads, 
right  to  left  as  viewed  facing  the  edge  connector.  The  signal  levels  at  this  point  in  the 
circuit  may  be  monitored  using  the  test  loops  adjacent  to  these  intermediate  amplifiers, 
aligned  with  the  longitudinal  axis  of  the  board. 

2.  Output  Buffer  Amplifiers 

For  convenience  in  recording  and  monitoring  requirements,  the  intermediate 
amplifier  outputs  are  also  connected  via  buffer  amplifiers  (U8A,  USB,  and  U8C  in  Figure 
Bl)  to  the  lower  row  of  BNCs  available  on  the  front  panel.  These  correspond  to  channel 
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numbers  1-3  (as  numbered  on  the  left  hand  edge  of  the  board)  from  left  to  right,  viewed 
facing  the  front  panel. 

C.  MODULATION  DEPTH  ADJUSTMENT 

The  modulation  depth  detected  in  each  fiber  has  been  observed  to  vary  slightly 
between  legs.  Adjustment  of  the  AC  portion  for  equal  signal  level  between  channels  also 
changes  the  DC  portion  of  the  signal.  A  set  of  summing  amplifiers  has  been  included  to 
remove  the  DC  portion  of  the  signal.  This  approach  an  alternative  to  high-pass  filtering  the 
input  signal  since  sub-radian  interferometric  signals  contain  DC.  This  approach  has  proven 
successful,  provided  that  the  ratio  of  modulation  depth  variations  between  channels 
remains  constant.  Changes  in  the  ratios  of  3x3-coupler  outputs  seems  to  cause  problems 
Investigation  of  the  demodulator  performance  is  required  to  address  interchannel 
variations. 

The  adjustment  to  remove  the  DC  portion  of  the  signal  is  made  via  the  potentiometers 
(tan  colored)  near  the  middle  of  the  board.  The  effect  of  this  adjustment  to  each  channel 
can  be  monitored  using  the  test  loops  under  the  differentiator  resistors  (mounted  on 
component  towers  in  the  middle  of  the  board).  The  signals  at  these  test  points  should  be 
adjusted  such  that  they  have  equal  AC  amplitudes  and  zero  DC  component.  Adjustment  to 
correct  discrepancies  should  be  performed  using  the  input  amplifier  potentiometers  (blue) 
first  and  then  the  modulation  depth  potentiometers  (tan). 

D.  DIFFERENTIATOR  ADJUSTMENT 

The  signal  presented  to  the  multiplier  sections  should  be  as  large  as  possible  without 
distortion  (10  V).  To  this  end  the  differentiator  transfer  function  at  the  highest  frequency 
of  interest  should  be  unity.  The  resistors  mounted  on  the  component  headers  permit 
adjustment  of  this  transfer  function.  In  this  implementation  40.2  kO  resistors  are  fitted. 
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corresponding  to  an  upper  frequency  of  interest  of  approximately  40  kHz,  as  computed  in 
the  following  formula. 


1 

2^RC 


_ 1 _ 

2;r(40.2knX100pF) 


39  6  kHz, 


(Bl) 


E.  PERFORMANCE 


Reference  9,  Chapter  XI,  details  the  performance  of  the  first  implementation  of  the 
symmetric  analog  demodulator  (SAD);  Reference  9  also  details  a  digital  implementation  of 
an  Asymmetric  Digital  Demodulator.  Following  the  analysis  of  Reference  9,  page  161,  the 
analog  interferometric  simulators  were  used  to  provide  inputs  to  the  SAD  V2.1.  This 
removes  the  added  noise  due  to  laser  phase  and  the  wandering  or  drifi  of  the  polarization 
of  the  optical  waves  in  the  interferometer,  the  major  cause  of  drift  in  amplitude  of  the 
interferometric  signal.  The  Scale  Factor  of  the  demodulator  is  defined  as  the  ratio  of  the 
demodulator  output  voltage  modulation  amplitude  to  the  phase  modulation  amplitude  of 
the  input  interferometric  signals. 

The  phase  noise  for  the  demodulator  at  5  kHz  was  measured  as  6.33  prad/VHz  The 
average  phase  noise  spectral  density  of  the  three  analog  interferometric  simulators  was 
measured  at  1  kHz  to  be  0.48±0.005  prad/>/Hz  [Reference  15,  book  1,  page  19]  and 
although  not  measured  at  5  kHz  is  still  far  enough  below  the  phase  noise  floor  of  the 
demodulator  that  it  would  not  be  of  the  order  of  the  demodulator  noise  at  5  kHz.  Thus  the 
demodulator  noise  measurement  can  be  mostly  attributed  to  the  demodulator  itself  The 
dynamic  range  at  5  kHz  was  determined  for  maximum  signal  (4%  THD)  to  be; 

=  7.4x10* 

6.33 /a-ad  (B2) 

Dynamic  Range  at  5  kHz  =  201og(7.4x  10*)  =  137  dB. 
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APPENDIX  C.  DEMONSTRATION  WITH  ANALOG  ELECTRONICS 


Prior  to  an  optical  system  demonstration  a  single  channel  demonstration  of  the 
multiplexing  technique  was  attempted  for  proof  of  concept  purposes.  Since  a 
demonstration  with  only  a  single  interferometric  output  from  the  3x3-coupler  was 
attempted,  the  demodulation  of  the  interferometric  signals  could  not  be  carried  out.  The 
essence  of  this  experiment  was  that  if  the  ratios  of  the  demultiplexed  components  of  the 
interferometric  signal  and  the  original  signal  remained  constant  throughout  the  various 
steps  of  frequency  translation,  then  one  could  conclude  that  the  technique  was  feasible 
That  is  to  say,  if  the  fundamental  and  significant  higher  order  harmonics,  which  are  greater 
than  the  noise  floor,  were  in  the  same  ratios  as  the  original  interferometric  signal, 
demodulation  would  produce  the  same  output  phase  modulation  for  both  the  baseband 
and  the  demultiplexed  signals.  Further,  the  demonstration  did  not  quantify  the  effects  of 
the  loss  of  higher  order  harmonics  (as  a  result  of  the  multiplexing  scheme).  The  effect 
would  be  extremely  small  but  could  be  quantified  in  terms  of  the  amount  of  energy  lost 
between  the  baseband  signal  and  the  demultiplexed  signal.  This  is  easily  computed  from 
the  fact  that  the  energy  contained  in  the  interferometric  signal  that  can  be  demodulated  is 
the  sum  of  the  amplitude  squared  values  of  all  of  the  components  (Parseval's  theorem).  By 
calculating  the  energy  before  and  after  multiplexing,  and  the  energy  lost,  a  reasonable 
estimate  of  the  effect  of  the  distortion  fi'om  the  demodulated  signal  can  be  made.  It  will  be 
shown  later  that  this  effect  is  small.  For  low  level  signals  however,  it  provides  a  lower 
limit  on  the  level  of  signal  which  can  be  demodulated,  i.e.,  the  minimum  detectable  signal. 
This  analysis  assumes  that  dispersion  of  the  higher  order  harmonics  is  minimal,  i.e.,  the 
higher  order  harmonics  are  integral  multiples  of  the  fundamental. 
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Another  aspect  of  the  demultiplexing  technique  not  covered  in  this  demonstration  is 
the  effect  of  filtering  on  the  interferometric  signal  prior  to  demodulation  This  was  not 
considered  since  filter  theory  is  well  understood  and  documented,  and  the  only 
requirements  to  ensure  successful  demodulation  are  constant  gain  and  phase  over  the  pass 
band  of  interest.  This  passband  of  interest  is  the  interferometric  signal  significant 
bandwidth/,^  discussed  earlier  (section  II  B  I  .). 

The  equipment  setup  for  the  demonstration  is  shown  at  Figure  Cl . 
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Figure  Cl.  Equipment  Setup  for  Aiinlog  Electronics  Dciuoiistrntion. 
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The  setup  has  three  separate  function  generators  to  simulate  the  acoustic  signal  from 
the  sensor  and  wavelength  and  intensity  modulation  from  the  laser  source  The  wavelength 
modulation  is  phase  and  frequency  locked  to  the  intensity  modulation  of  the  source  via 
appropriate  operation  of  the  triggers  of  the  respective  function  generators.  Another 
function  generator  is  utilized  as  the  demultiplexer  or  mixer  by  appropriate  use  the  AM 
input.  The  outputs  were  monitored  by  a  HP3562A  spectrum  analyzer. 

A.  ANALOG  INTERFEROMETRIC  SIMULATOR 

The  analog  interferometric  simulator  is  similar  to  that  described  at  Reference  9, 
Appendix  C.  A  modification  was  conducted  so  that  the  input  previously  labeled  SENSOR 
2  is  now  labeled  STATIC  PHASE  0/P  and  provides  an  easier  way  to  measure  the  static 
phase  drift  angle  Previously  this  was  determined  using  Lissajous  figures  [Reference 
20,  pages  30-39]  and  making  appropriate  measurements  [Reference  9,  Appendix  C, 
section  B,  pages  301-2].  The  drift  angle  term  in  the  simulator  is  not  time- varying  as  it 
would  be  in  an  actual  system  but  a  constant  set  on  the  front  panel  via  the  static  phase 
adjustment. 

The  modified  simulator  block  diagram  is  shown  in  Figure  C2  while  the  circuit  is 
shown  in  Figure  C3.  The  static  phase  measurement  calculations  will  be  presented  by 
isolating  that  part  of  the  circuit  involved,  see  Figure  C4 
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I'igiirc  C2.  Modified  Analog  liitciTcronietnc  Signal  Sininlator  Block  Diagram. 


Figure  C3.  Modified  Analog  liiiviferometric  Signal  Simulator  Circuit  Diagram. 
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Tlie  Analog  Devices  AD639  Universal  Trigonometric  Function  Converter  is 
conllgurcd  to  produce  the  sine  of  its  input.  A  one-volt  input  is  equivalent  to  an  input 
argument  of  50°  Thus  the  gain  G  of  the  circuit  of  Figure  C4  is: 


G 


2 


2.85. 


(Cl) 


Thus  the  static  phase  adjustment  input  equates  to; 

2.85  x  50°/Voll  =  l42°/Volt.  (C.2) 


so  that,  with  one  simulator  adjusted  to  zero  volts  static  phase  output  and  the  other  two 
static  phase  outputs  adjusted  to  ±0.845  V,  we  obtain  our  three  interferometric  signals 
symmetrically  phase  shifted  120°  relative  to  each  other. 


STATIC  PHASE 
MEASUREMENT 


Figure  C4.  Circuit  for  Static  Phase  Measurement  Calcuiatioiis. 
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B.  RESULTS 

The  first  test  verified  that  the  modulated  components  were  in  the  same  ratio  as  the 
baseband  components.  The  input  labeled  PZT MOD  in  Figures  Cl  and  C2  was  grounded 
to  simulate  the  absence  of  wavelength  modulation.  The  static  phase  was  set  at  120° 
(0.845V  static  phase  output)  to  ensure  that  both  even  and  odd  harmonics  would  be 
present  in  the  interferometric  signal.  A  1  kHz  sine  wave  was  applied  to  the  SENSOR  1 
input  to  generate  an  interferometric  type  signal  The  amplitude  was  varied  and  the  ratios 
of  baseband  components  and  modulated  components  were  observed  to  'track'  each  other 
A  plot  of  a4.2Vpp  amplitude  =  2. 1  radians  from  simulator  scale  factor  of  1  rad/V)  is 
shown  at  Figure  C5  with  the  upper  waveform  showing  the  baseband  'interferometric 
signal'  prior  to  multiplexing,  and  the  lower  signal  showing  the  amplitude  modulated  and 
demultiplexed  version. 
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Figme  C5.  Speciruiti  of  baseband  simulated  interrerometric  signal  (top)  and 
speciriiiii  of  demultiplexed  version  of  same  inlerferometnc  signal 
(bottom). 


It  can  be  seen  that  signal  and  demultiplexed  components  are  indeed  in  the  same  ratios 
even  though  the  amplitudes  are  different . 

Figure  C6  shows  the  Lissajous  pattern  for  the  demultiplexed  output  (x  axis)  and  the 
simulated  interferometric  signal  output  (y  axis)  with  no  wavelength  modulation  (PZT 
MOD  input  grounded).  This  also  verifies  that  the  signals  are  in  the  same  phase  and  of  the 
same  shape,  i.e  ,  the  component  ratios  are  the  same.  The  closed  ellipse  or  nearly  straight 
line  offset  at  45°  indicates  a  zero  phase  difference.  The  total  harmonic  distortion  THD  is 
the  same  for  both  waveforms.  The  THD  together  with  Figure  C5  also  provides  a  crude 
measure  of  the  ratios  of  components. 

The  HP3562A  Spectrum  Analyzer  calculates  the  THD  by  calculating  the  ratio  of 
energy  in  the  fundamental  component  to  the  sum  of  the  other  components  Thus  while  it  is 
possible  to  have  the  components  in  different  ratios  and  the  same  THD,  it  is  not  possible  to 
have  them  in  different  ratios  and  the  same  THD. 
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Wavelength  modulation  of  the  source  is  simulated  by  a  variable  amplitude  signal 
locked  in  phase  and  frequency  to  the  modulation  frequency  connected  to  the  PZT  MOD 
input.  This  wavelength  modulation  was  not  included  in  Equation  13  (where  it  was 
presumed  negligible).  Equation  13  is  rewritten  below  with  the  wavelength  modulation 
included. 


u  =  A  +  Bcos(^t)  +  ^j^(t)+^, +^„sinn>„t).  (C.3) 


Hence  represents  the  amount  of  wavelength  modulation.  Figure  C7  shows  the 
Lissajous  pattern  for  the  case  where  (j>,„  =  0.42  radians  for  an  'acoustic'  signal  of  2  12 
radians  amplitude  at  1  kHz  signal  and  a  modulation  frequency  yj,,  =  20  kHz  This  clearly 
shows  that  the  two  signals  are  no  longer  in  phase  (open  ellipse)  or  the  same  shape  (seen  as 
the  jagged  shape  of  the  outside  of  the  ellipse).  This  is  further  highlighted  by  Figure  C8 
where  the  THD  ratio  and  component  ratio  are  no  longer  the  same 

This  change  highlights  the  difference  between  this  multiplexing  scheme  and  the 
presently  popular  wavelength  modulated  frequency  division  multiplexing  (WM-FDM) 
WM-FDM  is  more  commonly  referred  to  as  phase  generated  carrier  (PGC)  [Reference  5, 
section  11.4.2.2,  pages  360-2  and  Reference  1,  volume  II,  section  14.6.5,  pages  564-5]. 
While  WM-FDM  seeks  to  maximize  the  wavelength  modulation  and  to  reduce  the 
intensity  modulation,  our  technique  AM-FDM  operates  optimally  with  no  wavelength 
modulation. 

The  analog  electronics  demonstration  showed  that  the  multiplexing  technique  will 
work  provided  that  wavelength  modulation  is  not  significant.  The  extent  to  which  the 
proposed  multiplexing  technique  will  work  with  increasing  or  decreasing  wavelength 
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luodiilalion  has  not  been  studied.  When  the  amount  of  wavelength  modulation.  is 
approximately  3%  of  the  signal  amplitude.  the  Lissajous  pattern  starts  to  change  from  a 
line  to  form  an  open  ellipse. 


Figure  C7.  Lissajous  pattern  for  same  signal  as  Figures  C5  and  C6  with  30% 
waveleuglli  modulation  0.300  rad). 
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Figure  C8.  Spectriiiii  of  signals  of  Figure  C7  (with  wavelength  modulation); 

baseband  interreronietric  signal  (top)  and  spectrum  of  demultiplexed 
version  of  same  inteiTerometric  signal  (bottom). 
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APPENDIX  D.  MISCELLANEOUS 


This  appendix  contains  the  following; 

•  MATLAB™  routines  written  by  the  author  for  this  thesis. 

•  THB292  22  kHz  integrated  active  filter  response. 

•  Operating  curves  and  data  sheets  for  lasers  used  in  this  thesis. 

•  Benchtop  Interferometer  Figure. 

•  ORTEL™  information  sheet  on  device  linearity  specifications. 

•  Noisefloor  Figures  for  Analog  Interferometric  Simulator  and  Mixing  circuit 

The  MATLAB™  routines  are  provided  with  no  guarantee  of  results  obtained  i  e.,  the 
author  disclaims  any  errors  in  the  code  reproduced. 
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MATLAB™  routine  used  for  Chapter  II  figures. 


%  Thciis  -  Mixing  of  intoferomcthc  signal  simulation 
clear.clg; 

6 1 » lOOO^ihis  1 

62-600ailus2”pi;phkl2-0; 
finl >40000; 
fin3>60000; 
phia>pi; 

%  phid  is  the  drift  phase  term  for  the  difference  in  optical  path  length 

%  between  the  two  legs  of  the  uiterferomeler 
%iphid>input('£nter  phid  (in  radians)  *); 
t-0:(  10/(61  •32767));(10/&1 ); 

%l=0;(  10/(fs»32767));(  10/fs); 
m=0.75. 

1 1 «( 1  +(cos(phisl  •sin(2*pi*61  •t)  +phidl  ))).*(  1  +(m«sin(2»pi»lmi»t))); 
I2=(  1  +(cos(phis2*sin(2*pi*fs2»t)  +phid2))).*(  1  ^^m*sin(2*pi»fin2»l))); 

I3=(n+I2); 

I=I3.*sin(2*pi*fml*t); 

F=ffl(1.32768). 

f=0;  l/iiiax(t):32767*  l/max(t). 

(k=f./1000; 

subplot(2 1  IXplotlff  ( 1 : 1 200),20*logl  0(abs(F(  1 ;  1 200)))).grid 
xlabeK'FrequeiKy  (kH2)'),ylabelCMagnitude  (dB)') 
axis((0  fled  200)0  100)); 
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MATLAB™  routine  used  for  Chapter  III  figures. 


SThcai*  nlculatiom  for  MxN  array 
clear,ctg; 

N-inpitfCEiiter  #  of  icnaon  N  *); 
M'lnputCEnlcr  *  of  icmaci  M  7i 


fori-l:N 


I(iH(  I  - 1  1  -(iyN>+  I/Ny  J  .230268770812; 

c(i)=((  1  - 1 /N  VN  V((  1  - 1 /N  )•(  1 )» I /N  ); 
ptt(i)=l(i)»c(i); 

pkl(iH).003*I(i)*e(i); 

end 

c' 

fdB=  1 0*logl  0(anes(N,  1  )-c’) 

pkdb=  1 0*log  1 0(pk);pk  1  db*’ 1 0*log  1 0(pk  1 ), 

%:  I  ^flipud(c<  I  :lengUt(c>-l )'). 

•/«forj=2:lenglh(c>-I 

O/Kind 


for  i=l:M 


Iin(i)={(  1  -1  /M)*(  1/My  1 . 1 455 1 2941446; 

cni(i)=((l-l/M)/M)/((l-l/M)*(l-i/M)*l/M); 

•/.pk(i)=l(i)*c(i); 

•/ipkl(i)=0.003»l(i)*c(i); 

end 

cin' 

fdBin=  10*logl0(ones<M,  1  )<ni') 

%pkdb=  10*iogl0(pk);pk  I  db=  1 0’logl  0(pk  I ); 
c  1  =f]ipud(cm(  1  :length(cni)- 1 )'); 
ll(l)=plcl(l); 


for  i-2:lenglh(cm) 

Il(i)-(cl{i-l)Ml(i.l)); 

end 

IIdB°10*logl0(Il); 

figurefl) 

plol(I,'o7 

ylabel(Tractional  Power  Input  Pin/Ptotal') 
vlabelCSentor  #  N'Xgrid 


figuirv(2) 

plo«(c,V).yUbeK‘Coupl«'  RaJio) 

KlabeK'Coupkr  » (N-1  toul)').gn<i 

sxu((ON  0  0.61) 

fif{ure<3) 

plo<(fda'o'),yUbel(X^pta'  Raiio  (dB)') 
xIabeU'Coupler  tt  (N-1  U)Ul)’).grid 
axis((0  N  -5  01) 

figure(4) 

plot(pkdb,'o‘) 

yUbel('Po«**r  to  Sensor  N  Pn/PloUl  (dB)') 
xIabelCScnsor  *»  N‘),grid 

figure(}) 

plol(pkldb-2.04,'o'),ylabel('Pow«f  Sensor  N  Po“3niW  (dBin)’) 
xlabeK'Sensor  4  N'),grid 

rigure<6) 

plol(cl.'o'),grid 

ylabeICCoupliiig  Ralio'),xlabe)CCoupler  »  (M-1  toul)') 
figure(7) 

plot(l  I  dB-2.04,'o‘),gnd 

ylabel('Power  at  Sensor  M  (dB)  ),xlabel('Sensor «  M') 


MATLAB™  routine  used  for  figure  A2. 


%  IiHetfcronielhc  Signal  Simulation 
S 

H  for  A/B=  1  1(1)- 1  +ooi(  phis  (sin|oniegas  t]  +  phid)  J 
% 

clcar.clg; 

&-I000; 

%  phis  is  the  amplitude  of  the  signal  of  interest 

phis-input(^ter  phis  (in  radians)  *); 
Sphisl-inpulCEnter  phisl  (in  radians)  *); 
Sphis2-input('Enter  phis2  (in  radians) '); 


%  phid  is  the  drift  phase  term  for  the  difference  in  optical  path  length 
%  between  tlie  two  legs  of  the  interferometer 
phid-input('£nter  pliid  (in  radians)  ‘); 
phidl-input('Enter  phid  (in  radians)  '); 
phid2-input(‘Enter  phid  (in  radians)  ‘); 

t=0.(l0/(fs*32767)).(10/fs); 

1=  I  •t-(cos(pliis*sin(2*pi*fs*t)  -rphid)); 

11  =  1  +(cos(plus*sin(2*pi*fs»t)  -t-pltidl )), 

12= 1  +(cos(pliis*sin(2*pi*fs*t)  •>-phid2)); 

F=fll(1, 16384);  Fl=ffl(1 1.1 6384);  F2=ffl(12.16384); 
f-0:2/max(t);  16383*1  /ntax(l); 

subplot(221),plot(t(l:length(t)/1000:lengtl*'t)'5),I(l:lenglh(t)/1000:length(t)'5)),grid 
xlabel(Time  (sec)'),ylabel('lntensit>’‘) 
axis((0 1(70)0  100|); 

subplot(222),plot(l(1.400),20*logl0(abs(F(l:400)))),grid 
xlabel(Trequency  (Hz)'),ylabel(‘Magnitude  (dB)') 

$ubpIot(223),plot(t(l:length(t)/l000:lengtli(t)/5),n(l:length(t)/1000:length(t)'S)),grid 
xlabeKTime  (sec)'),ylabel(’lnlensity’) 
axis((0 1(70)0  100]); 

subplot(224).plot(f(  1 :200),20*logl  0(abs(F  1  ( 1 :200)))),grid 
xlabelCFrequency  (Hz)'Xylabel('Mapiitude  (dB)') 
pause 
cig; 

subplot(22I),plot(t(l:lenglli(tyiOOO:length(ty)),I2(t:length(iyiOOO:length(tV})),grid 

xlabelfTime  (secyXylabelCIntensity*) 

pause 

axis([0  f(70)0  100]); 

subplot(222Xplot(i(  1 :4000),20*logl  0(abs(F2(  1 :4000)))Xgrid 
xlabelCFrequency  (Hz)'),ylabelCMagnitude  (dB)') 
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I.ASER  DIODE  TEST  DATA 


Mo  d  e  I  No  .  : 
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Wa  va  I  a  n  g  t  H 
Ua  s  a  r  Type 
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Figure  D2.  Seaslar  Laser  #SP2676  Maiiiiractiirer's  lest  data. 
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UDT020  Power  0/P  (mW) 


Seastar  Laser  fSP2676 


Forward  Current  (mA) 


Figure  D3.  Seastar  Laser  #SP2676  measured  operating  curve. 
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LASER  DIODE  TEST  DATA 


T4odo  I  No  .  s 
Se  r  I  o  I  No . 
Wo  VO  I  o  ng  -t  h 
Laser  Type 


PT- AS0 
1  320 
B30  nm 

HITACHI  HLa312G 


Ra  tod  Po ! 
Oa  t  e  teste 
Tasted  by: 


1 d37  micro  W 


powen  OUTPUT  us.  roawAno  cunneuT 


"IN  CUnnSNT  vs.  OUTPUT  POWEn 


•  sn  inn  mn  ann  ann 

ronwAno  cunneuT  f«A,nc> 


4 


Figure  D4.  Seustnr  Lnscr  #SPI320  Maiiiirncliirer’s  lest  data. 


Power  (mW) 


Power  (mW) 


Drive  Current  (mA) 


Figure  D5.  Senslar  Laser  //SP1320  measured  operating  ciii've. 
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Figure  21 

Third  Order  tnlermodulalion 
Products  end  TNrd  Order 
Intercept  Llr>earHy  -  Two  Tone 
Test 
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prodiiTlS  will  be  mire  at  rniich 
below  the  earner  leivl  os  the 
carrier  is  below  the  input  third 
order  inlereepi  point  (TOI). 


Figure  D7,  ORTEL"^  device  linearity  speciricalions. 
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